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Abstract —  This  paper  describes  the  development 
and  verification  of  a  six  degree  of  freedom,  non-linear 
simulation  model  for  the  REMUS  AUV,  the  first  such 
model  for  this  platform.  In  this  model,  the  external 
forces  and  moments  resulting  from  hydrostatics,  hy¬ 
drodynamic  lift  and  drag,  added  mass,  and  the  con¬ 
trol  inputs  of  the  vehicle  propeller  and  fins  are  all 
defined  in  terms  of  vehicle  coefficients.  This  paper 
briefly  describes  the  derivation  of  these  coefficients. 
The  equations  determining  the  coefficients,  as  well  as 
those  describing  the  vehicle  rigid-body  dynamics,  are 
left  in  non-linear  form  to  better  simulate  the  inher¬ 
ently  non-linear  behavior  of  the  vehicle.  Simulation 
of  the  vehicle  motion  is  achieved  through  numeric  in¬ 
tegration  of  the  equations  of  motion.  The  simulator 
output  is  then  verified  against  vehicle  dynamics  data 
collected  in  experiments  performed  at  sea.  The  sim¬ 
ulator  is  shown  to  accurately  model  the  motion  of  the 
vehicle.  The  paper  concludes  with  recommendations 
for  future  model  validation  experiments. 
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I.  Introduction 

IMPROVING  the  performance  of  modular,  low- 
cost  autonomous  underwater  vehicles  (AUVs)  in 
such  applications  as  long-range  oceanographic  sur¬ 
vey,  autonomous  docking,  and  shallow-water  mine 
countermeasures  requires  improving  the  vehicles’ 
maneuvering  precision  and  battery  life.  These  goals 
can  be  achieved  through  the  improvement  of  the 
vehicle  control  system.  A  vehicle  dynamics  model 
based  on  a  combination  of  theory  and  empirical  data 
would  provide  an  efficient  platform  for  vehicle  con¬ 
trol  system  development,  and  an  alternative  to  the 
typical  trial-and-error  method  of  vehicle  control  sys¬ 
tem  field  tuning.  Furthermore,  a  good  vehicle  dy¬ 
namics  model  is  a  crucial  element  of  the  Kalman 
filter  at  the  heart  of  any  navigation  algorithm.  As 
there  exists  no  standard  procedure  for  vehicle  mod¬ 
eling  in  industry,  the  simulation  of  each  vehicle  sys¬ 
tem  represents  a  new  challenge. 

II.  The  REMUS  AUV 

Developed  by  von  Alt  and  associates  at  the 
Woods  Hole  Oceanographic  Institute’s  Oceano¬ 
graphic  Systems  Laboratory  (WHOI  OSL),  RE- 

This  work  was  supported  by  the  Office  of  Naval  Research 
under  grant  N0001 4-97- 1-0787. 


MUS  (Remote  Environmental  Monitoring  Unit)  is 
a  low-cost,  man-portable  AUV  design  with  approx¬ 
imately  1000  hours  of  water  time  over  hundreds  of 
missions  on  15  vehicles,  with  applications  in  au¬ 
tonomous  docking,  long-range  oceanographic  sur¬ 
vey,  and  shallow- water  mine  reconnaissance  [14], 
[15],  [16], 

The  hull  shape  of  the  REMUS  vehicle  is  based 
on  the  Myring  hull  profile  equations  [11],  which  de¬ 
scribe  a  body  contour  with  minimal  drag  coefficient 
for  a  given  fineness  ratio  (body  length/maximum 
diameter).  For  reference,  Myring  [11,  p.  189]  as¬ 
sumes  a  total  body  length  of  100  units,  and  classi¬ 
fies  body  types  by  a  code  of  the  form  a/b/n/0/^d , 
where  6  is  given  in  radians.  REMUS  is  based  on  the 
Myring  B  hull  contour,  which  is  given  by  the  code 
15/55/1.25/0.4363/5.  See  Figure  1  for  a  plot  of 
the  REMUS  profile,  and  Figure  2  for  a  picture  of 
the  vehicle. 


Fig.  1:  REMUS  Profile  (XZ-plane) 


Fig.  2:  The  REMUS  AUV 


III.  Vehicle  Model  Development 

The  standard  submarine  equations  of  motion — 
developed  by  Gertler  and  Hagen  [7]  and  revised  by 
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Humphreys  [10]  and  Feldman  [4] — offer  a  general 
framework  for  the  development  of  the  vehicle  equa¬ 
tions  of  motion.  In  these  equations  of  motion,  ex¬ 
ternal  forces  and  moments 

^  ^  Fc,t  =  -^hydrostatic  H”  ^lift  HH  ^drag  T  -Fcontrol 

are  described  in  terms  of  vehicle  coefficients.  For 
example,  vehicle  axial  drag: 

Fd  =  -  Q/rcdJ4/)  u  |«|  =  Xu\u\u  |w| 

where  the  coefficient 

„  9Fd  _  1  4 

XuW~  d(u\u\)~  2pCdAf 

The  vehicle  coefficients  were  derived  as  follows: 

•  Axial  Drag  Author  tow  tank  experiments  [1]. 

•  Added  Mass  Hydrodynamic  theory  from  New¬ 
man  [12]  and  empirical  formulae  from  Blevins  [2]. 

•  Vehicle  Crossflow  Drag  Hydrodynamic  theory 
from  Newman  [12]  and  empirically-derived  formu¬ 
lae  from  Hoerner  [8]. 

•  Vehicle  Body  Lift  Empirically-derived  formulae 
from  Bottaccini  [3] ,  Fidler  [5]  and  Hoerner  [9] . 

•  Fin  Forces  Empirically-derived  formulae  from  Ho¬ 
erner  [8],  [9]  and  Whicker  and  Fehlner  [17]. 

The  author  uses  a  very  simple  model  for  the  RE¬ 
MUS  propulsion  system,  which  treats  the  propeller 
as  a  source  of  constant  thrust  and  torque.  The 
values  for  these  coefficients  are  based  on  a  vehicle 
forward  speed  of  three  knots  (1.54  m/s),  and  were 
derived  from  vehicle  design-stage  propeller  bench 
tests  conducted  by  Ben  Allen  at  WHOI  OSL,  and 
from  experiments  conducted  by  the  author  at  sea. 
This  simple  model  is  acceptable  for  small  amplitude 
perturbations  about  the  vehicle  steady  state.  In 
the  future,  the  author  would  like  to  replace  this 
with  a  more  sophisticated  propeller  model,  such 
as  developed  by  Yoerger  and  Slotine  [18],  or  with 
experimentally-derived  values. 

IV.  Modeling  Assumptions 

The  following  assumptions  were  made  in  the  de¬ 
velopment  of  the  vehicle  model: 

•  The  vehicle  is  a  rigid  body  of  constant  mass.  In 
other  words,  the  vehicle  mass  and  mass  distribution 
do  not  change  during  operation. 

•  Control  surface  assumptions.  We  assume  that  the 
control  fins  do  not  stall  regardless  of  angle  of  at¬ 
tack.  We  also  assume  an  instantaneous  fin  response, 
meaning  that  that  vehicle  actuator  time  response  is 
small  in  comparison  with  the  vehicle  attitude  time 
response. 


•  The  vehicle  is  deeply  submerged  in  a  homoge¬ 
neous,  unbounded  fluid.  In  other  words,  the  vehicle 
is  located  far  from  free  surface  (no  surface  effects, 
i.e.  no  sea  wave  or  vehicle  wave-making  loads),  walls 
and  bottom. 

•  The  vehicle  does  not  experience  memory  effects. 
The  simulator  neglects  the  effects  of  the  vehicle 
passing  through  its  own  wake. 

V.  Equations  of  Motion 

The  vehicle  equations  consist  of  the  following  el¬ 
ements: 

•  Kinematics:  the  geometric  aspects  of  motion 

•  Rigid-body  Dynamics:  the  vehicle  inertia  matrix 

•  Mechanics:  forces  and  moments  causing  motion 
For  a  treatment  of  the  generic  vehicle  kinematic 
equations,  see  Fossen  [6]. 

Combining  the  equations  for  the  vehicle  rigid- 
body  dynamics  with  the  equations  for  the  forces  and 
moments  on  the  vehicle,  we  arrive  at  the  combined 
nonlinear  equations  of  motion  for  the  REMUS  vehi¬ 
cle  in  six  degrees  of  freedom.  These  equations  follow 
the  SNAME  convention  for  the  assignment  of  the 
body-fixed  vehicle  coordinate  system. 

Surge,  or  translation  along  the  vehicle  x-axis: 

(m  -  Xu)u  +  mzgq  -  mygr  —  Xhs  +  Xu\u\u\u\ 

+  (Xwq  -  m)wq  4-  ( Xqq  +  mxg)q 2  ^ 

4*  ( Xvr  4-  m)vr  +  ( Xrr  +  mxg)r 2 
—  mygpq  —  mzgpr  +  Xpr0p 

Sway,  or  translation  along  the  vehicle  y-axis: 

(m  —  Yij)v  —  mzgp  4-  (mxg  -  )r  =  Yhs 
+  Y„Mt,  M  +  YrWr  |r|  +  mygr2 
+  (Yur  —  m)ur  4-  ( Ywp  +  m)wp  (2) 

4-  ( Ypq  —  mxg)pq  4-  Yuvuv  4-  mygp 2 
4-  mzgqr  4-  YuuSr  u2Sr 

Heave,  or  translation  along  the  vehicle  z-axis: 

(m  —  Z^)w  4-  mygp  -  (mxg  4-  Zg)q  =  ZHs 
+  Zw\w\w\w\  +  Zq\q\q\q\ 

4-  ( Zuq  +  m)uq  4-  ( Zvp  -  m)vp  (3) 

4-  ( Zrp  —  mxg)rp  4-  Zuwuw  4-  mzg(p2  4 -  q2) 

-  mygrq  4-  Zuus3u2S3 

Roll,  or  rotation  about  the  vehicle  x-axis: 

-  mzgv  4-  mygw  4-  (/**  -  Kp)p  =  KHs 

+  ^p|P|P|p|  —  (/«  “  Iyy)Qr  (4) 

4-  m(uq  —  vp)  —  mzg (wp  —  ur )  4-  -Kprop 
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Pitch,  or  rotation  about  the  vehicle  y-axis: 

mZgU  —  ( mxg  +  Mw)w  +  {Iyy  —  Mq)q  =  MhS 
+  MW\W\W  |w|  +  Mq\q\q  | <7 1 

4  (Muq  —  mxg)uq  4  ( Mvp  +  mxg)vp  (5) 

4  [Mrp  -  (Ixx  -  Izz) ]  rp  +  mzg(vr  -  wq) 

4  Muwuw  +  Muussu2Ss 

Yaw,  or  rotation  about  the  vehicle  z-axis: 

-  mygu  +  (mxg,  -  Nv)v  +  (?**  ~  N+)r  = 

+  ^«|t7|v|v|  +  Nr|r|r|r| 

+  (Nur  —  mxg)ur  4  ( Nwp  +  mxg)wp  (6) 

+  \Npq  -  (Iyy  ~  /**)]  P9  -  myg  (vr  -  w q) 

4-  Nuvuv  4  NuuSru2Sr 

Note  that  the  vehicle  cross-products  of  inertia, 
J^,  /**  and  Jyz,  are  assumed  to  be  small  and  are 
neglected  in  the  above  equations  of  motion.  Simi¬ 
larly,  the  equations  do  not  include  zero- valued  coef¬ 
ficients. 

VI.  Numerical  Integration 


Fig.  3:  The  author  (left)  and  Mike  Purcell  from 
WHOI  OSL,  running  vehicle  experiments  at  the 
Rutgers  Marine  Field  Station  in  Tuckerton,  NJ 
[Photo  courtesy  of  Nuno  Cruz,  Porto  University ] 


We  combine  the  above  equations: 

xn+i  =  xn  4-  —  (ki  -I-  2fc2  +  2&3  4-  k±)  (12) 


The  nonlinear  differential  equations  defining  the 
vehicle  accelerations  and  the  kinematic  equations 
give  us  the  vehicle  accelerations  in  the  different  ref¬ 
erence  frames.  Given  the  complex  and  highly  non¬ 
linear  nature  of  these  equations,  we  will  use  numer¬ 
ical  integration  to  solve  for  the  vehicle  speed,  posi¬ 
tion,  and  attitude  in  time. 

Consider  that  at  each  time  step,  we  can  express 
the  vehicle  equations  of  motion  as  follows: 

xn  =  f  (xn,Ttn)  (7) 

where  xn  is  the  vehicle  state  vector: 

xn  —  [uvwpqrxyz<j)6  ip]T  (8) 


to  yield  the  new  vehicle  state  at  each  time  step. 

VII.  Field  Experiments 

In  order  to  verify  the  accuracy  of  the  vehicle 
model,  the  author  conducted  a  series  of  experiments 
at  sea  measuring  the  response  of  the  vehicle  to  step 
changes  in  rudder  and  stern  plane  angle.  These  ex¬ 
periments  were  conducted  at  both  the  Woods  Hole 
Oceanographic  Institution  and  at  the  Rutgers  Uni¬ 
versity  Marine  Field  Station  in  Tuckerton,  New  Jer¬ 
sey  with  the  assistance  of  the  WHOI  OSL  staff. 

In  each  experiment  at  sea  we  measured  the  vehicle 
depth  and  attitude,  represented  in  the  vehicle  model 
by  the  following,  globally- referenced  vehicle  states: 


and  un  is  the  input  vector: 

y>n  —  [  &s  &r  Ypr0p  Kprop  ]  (9) 

We  will  use  the  Runge-Kutta  method  of  numerical 
integration,  for  which  we  first  calculate  the  follow¬ 
ing: 


fci  —  xn  4~  f  (xn,  un ) 


k2  =  f(x+  Yfel>“n+|) 

k3  =  f  (x  +  Yk2’un+h) 


&4  =  /  (®  +  Atks,un+1) 
where  the  interpolated  input  vector 

wn+£  =  2  ^Un  Wn+1) 


(10) 


(11) 


x  =  [  z  <t>  0  ]T  (13) 

In  these  experiments  we  also  recorded  the  vehicle 
fin  angles,  represented  in  the  vehicle  model  by  the 
following  vehicle-referenced  control  inputs: 

Un  =  [  6r  f  (14) 

Note  that  for  all  of  the  field  tests  described  in  this 
section,  the  vehicle  propeller  was  not  used  as  a  con¬ 
trol  input,  but  was  instead  kept  at  a  constant  1500 
RPM.  As  propeller  thrust  and  torque  were  difficult 
to  estimate  for  different  propeller  RPMs,  sticking 
to  a  constant  value  allowed  us  to  remove  a  source  of 
uncertainty  from  the  vehicle  model  comparison. 

In  order  to  measure  the  vehicle  response  to  step 
changes  in  fin  angle,  the  vehicle  was  given  the  fol¬ 
lowing  commands: 
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•  Timer  to  desired  depth  For  ” pitch  up”,  the  vehi¬ 
cle  was  commanded  to  six  meters  depth,  to  avoid 
breaking  the  surface.  For  "pitch  down”  commands, 
the  vehicle  was  commanded  to  2  meters  depth. 

•  Step  change  in  fin  angle  Upon  achieving  depth, 
the  vehicle  was  commanded  to  hold  a  certain  fin 
angle  for  two  seconds  in  the  case  of  vertical  plane 
response,  or  longer  for  horizontal  plane  response. 
The  fin  angle  duration  of  two  seconds  was  chosen  as 
a  result  of  the  experimental  run  shown  in  Figure  4. 
In  the  depth  plot,  right  around  the  seven-second 
mark,  you  can  see  that  the  vehicle  ran  into  and 
bounced  off  the  bottom.  Given  the  unpredictable 
vehicle  open-loop  response,  the  author  thought  it 
wise  to  use  short  periods. 


Fig.  4:  REMUS  Mission  Data:  Vehicle  bounces  off 
the  bottom  [d980729a,  Obj.  6) 

From  these  vehicle  experiments,  we  get  measure¬ 
ments  for  the  vehicle  response  to  step  changes  in 
rudder  and  stern  plane  angle.  It  is  important  to 
note  that  during  straight  and  level  flight,  the  vehi¬ 
cle  operates  at  a  roll  offset  of  negative  five  degrees 
(<j>  =  —5)  due  to  the  propeller  torque.  As  a  result, 
we  never  get  pure  vertical-  or  horizontal-plane  mo¬ 
tion.  That  said,  the  vehicle  roll  is  small  enough  that 
we  are  still  able  identify  the  vehicle  behavior  in  pitch 
and  yaw. 

See  Figure  5  for  REMUS  motion  while  operat¬ 


ing  under  closed-loop  control,  for  comparison  with 
the  open-loop,  step  response  data.  In  the  example 
shown,  the  vehicle  was  commanded  to  maintain  a 
depth  of  two  meters. 


Fin  Angles  vs  .  Ttne 


Vehicle  Depth  vs.  Time 


:::::: 

*___  :  1  — 

,  /  i / 

Vehicle  Yaw  vs.  Time 


Fig.  5:  REMUS  Mission  Data:  Vehicle  under  closed- 
loop  control.  [<1990727,  Obj.  4} 


VIII.  Model  Comparison 

The  model  was  given  initial  conditions  and  fin 
inputs  to  match  the  experimental  data.  The  follow¬ 
ing  uncertainties  affected  the  accuracy  of  the  model 
comparison: 

•  Vehicle  Initial  Conditions  The  greatest  uncer¬ 
tainty  was  the  vehicle  state  at  the  start  of  each  ex¬ 
perimental  objective.  We  were  unable  to  measure 
currents,  wave  effects,  and  non- axial  vehicle  veloci¬ 
ties,  which  would  have  all  affected  the  vehicle  mo¬ 
tion  during  open-loop  maneuvers. 

•  Control  Fin  Alignment  Although  the  alignment 
of  the  vehicle  fins  was  checked  before  each  exper¬ 
imental  mission,  it  was  difficult  to  keep  the  vehi¬ 
cle  control  fins  from  getting  knocked  during  vehicle 
transportation  and  launch.  This  could  have  resulted 
in  fin  misalignments  as  great  as  five  degrees. 

•  Attitude  Sensor  Dynamics  The  vehicle  attitude 
sensor  was  sensitive  to  coupling  due  to  vehicle  ac¬ 
celerations.  Although  most  likely  a  small  effect,  the 
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author  did  not  have  the  opportunity  to  characterize 
these  sensor  dynamics. 

Recommended  improvements  to  the  experimental 
method  and  instrumentation  are  discussed  in  the 
conclusion. 

Figures  6  and  7  show  the  vehicle  response  to  step 
changes  in  rudder  angle.  In  Figure  6,  the  vehicle 
was  given  zero  fin  inputs  for  ten  seconds,  then  four 
degrees  of  positive  rudder  for  25  seconds,  then  four 
degrees  of  negative  rudder  for  30  seconds. 

Figure  7  shows  a  direct  comparison  of  the  experi¬ 
ment  and  simulator  data  for  a  period  of  steady  yaw 
rate  (i.e.  no  initial  transients).  The  simulated  ve¬ 
hicle  yaw  rate  is  shown  to  be  a  very  close  match  to 
the  experiment.  Discrepancies  between  the  vehicle 
depth  rates,  and  vehicle  pitch  and  roll  angles  likely 
have  to  do  with  differences  in  the  simulator  initial 
conditions. 

Figures  8  and  9  show  the  vehicle  response  to  a 
step  change  in  pitch  fin  angle.  The  vehicle  was  given 
zero  fin  inputs  for  two  seconds,  then  four  degrees  of 
negative  pitch  fin  for  two  seconds.  Figure  9  shows 
a  vehicle  depth  change  of  roughly  0.5  meters  and  a 
pitch  change  of  twenty  degrees,  which  both  compare 
well  with  the  experimental  data.. 

IX.  Experimental  Recommendations 

The  vehicle  dynamics  data  collected  in  the  exper¬ 
iments  were  limited  in  the  number  of  vehicle  states 
recorded,  and  the  accuracy  of  those  measurements. 
This  made  it  difficult  to  judge  the  validity  of  the 
model  comparison. 

For  future  experiments,  the  author  has  aug¬ 
mented  the  standard  REMUS  sensors  with  an  in¬ 
ertial  measurement  unit:  the  Crossbow  DMU- 
AHRS  (Dynamic  Measurement  Unit — Attitude  and 
Heading  Reference  System).  This  instrument  out¬ 
puts  magnetic  orientation,  accelerations  and  angu¬ 
lar  rates  on  three  axes.  This  will  significantly  im¬ 
prove  our  ability  to  measure  the  vehicle  initial  con¬ 
ditions  in  particular,  and  the  vehicle  motion  in  gen¬ 
eral. 

Second,  even  with  improved  vehicle  instrumen¬ 
tation  it  is  still  important  to  understand  the  vehi¬ 
cle  steady-state  conditions.  AUVs  like  REMUS  can 
be  unstable  when  operating  without  control.  Us¬ 
ing  the  inertial  measurement  unit,  it  will  be  possi¬ 
ble  to  identify  the  propeller  RPM  and  fixed  fin  an¬ 
gles  which  result  in  straight  and  level  vehicle  flight. 
These  settings  will  then  be  used  at  the  start  of  every 
experimental  run. 


X:  Forward  Displacement  JmJ 


Fig.  6:  REMUS  Simulator  Data:  Vehicle  trajectory 
for  vehicle  response  in  the  horizontal  plane 


Fig.  7:  REMUS  Simulator  Data:  Comparison  plots 
for  vehicle  response  in  the  yaw  plane 
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X.  Conclusion 


Y2  PLOT 


SIMULATOR  XYZ  PM  -CU 


Y :  Lateral  Displacement  frr] 


The  vehicle  dynamics  model  has  been  shown 
to  accurately  simulate  the  motion  of  the  REMUS 
AUV,  to  the  limit  of  the  existing  vehicle  sensors. 
For  a  detailed  description  of  the  model  development 
and  experimental  results,  please  refer  to  the  author’s 
MIT  Master’s  thesis  [13]. 

References 


XZ  PLOT 


2 1— : - = - ; - = . . 

0  1  2  3  4  5 

X:  Forward  Displacement  [m] 


XY  PLOT 


0  1  2  3  4  5 


X:  Forward  Displacement  [m] 


Fig.  8:  REMUS  Simulator  Data:  Vehicle  trajectory 
for  vehicle  pitching  up 


Time  (seconds) 


Fig.  9:  REMUS  Simulator  Data:  Comparison  plots 
for  vehicle  pitching  up 


[1]  B.  Allen,  W.  Vorus,  and  T.  Prestero.  Propulsion  system 
performance  enhancements  on  REMUS  AUVs.  In  Pro¬ 
ceedings  MTS/IEEE  Oceans  2000 ,  Providence,  Rhode 
Island,  September  2000.  2 

[2]  Robert  D.  Blevins.  Formulas  for  Natural  Frequency  and 
Mode  Shape .  Kreiger  Publishing,  Florida,  1979.  2 

[3]  M.  R.  Bottaccini.  The  stability  coefficients  of  standard 
torpedoes.  NAVORD  Report  3346,  U.S.  Naval  Ordnance 
Test  Station,  China  Lake,  CA,  1954.  2 

[4]  J.  Feldman.  Revised  standard  submarine  equations  of 
motion.  Report  DTNSRDC/SPD-0393-09,  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center, 
Bethesda,  MD,  June  1979.  2 

[5]  John  E.  Fidler  and  Charles  A.  Smith.  Methods  for  pre¬ 
dicting  submersible  hydrodynamic  characteristics.  Re¬ 
port  NCSC  TM-238-78,  Naval  Coastal  Systems  Labora- 
tory,  Panama  City,  FL,  1978.  2 

[6]  Thor  I.  Fossen.  Guidance  and  Control  of  Ocean  Vehi¬ 
cles .  John  Wiley  &  Sons,  New  York,  1994.  2 

[7]  M.  Gertler  and  G.  Hagen.  Standard  equations  of  mo¬ 
tion  for  submarine  simulation.  Report  DTNSRDC  2510, 
David  W.  Taylor  Naval  Ship  Research  and  Development 
Center,  Bethesda,  MD,  June  1967.  1 

[8]  Sighard  F.  Hoerner.  Fluid  Dynamic  Drag .  Published  by 
author,  1965.  2 

[9]  Sighard  F.  Hoerner  and  Henry  V.  Borst.  Fluid  Dynamic 
Lift.  Published  by  author,  second  edition,  1985.  2 

[10]  D.  E.  Humphreys.  Development  of  the  equations  of  mo¬ 
tion  and  transfer  functions  for  underwater  vehicles.  Re¬ 
port  NCSL  287-76,  Naval  Coastal  Systems  Laboratory, 
Panama  City,  FL,  July  1976.  2 

[11]  D.  F.  Myring.  A  theoretical  study  of  body  drag  in 
subcritical  axisymmetric  flow.  Aeronautical  Quarterly , 
27(3):  186-94,  August  1976.  1 

[12]  J.  N.  Newman.  Marine  Hydrodynamics.  MIT  Press, 
Massachusetts,  1977.  2 

[13]  Timothy  J.  Prestero.  Verification  of  a  six-degree  of  free¬ 
dom  simulation  model  for  the  REMUS  AUV.  Master’s 
thesis,  Massachusetts  Institute  of  Technology/Woods 
Hole  Oceanographic  Institution,  Departments  of  Ocean 
and  Mechanical  Engineering,  2001.  6 

[14]  M.  Purcell,  C.  von  Alt,  B.  Allen,  T.  Austin,  N.  Forrester, 
R.  Goldsborough,  and  R.  Stokey.  New  capabilities  of 
the  REMUS  autonomous  underwater  vehicle.  In  Pro¬ 
ceedings  MTS/IEEE  Oceans  2000 ,  Providence,  Rhode 
Island,  September  2000.  1 

[15]  C.  von  Alt,  B.  Allen,  T.  Austin,  and  R.  Stokey.  Re¬ 
mote  environmental  monitoring  units.  In  Proceedings 
MTS/IEEE  Oceans  1994 ,  Cambridge,  MA,  1994.  1 

[16]  C.  von  Alt  and  J.F.  Grassle.  LEO-15:  An  unmanned 
long  term  environmental  observatory.  In  Proceedings 
MTS/IEEE  Oceans  1992 ,  Newport,  RI,  1992.  1 

[17]  L.  F.  Whicker  and  L.  F.  Fehlner.  Free-stream  charac¬ 
teristics  of  a  family  of  low-aspect  ratio  control  surfaces. 
Technical  Report  933,  David  Taylor  Model  Basin,  1958. 
NC.  2 

[18]  D.  R.  Yoerger,  J.G.  Cooke,  and  J.-J.  E.  Slotine.  The 
influence  of  thruster  dynamics  on  underwater  vehicle 
behavior  and  their  incorporation  into  control  system  de¬ 
sign.  IEEE  Journal  of  Oceanic  Engineering ,  15:167-178, 
July  1990.  2 


455 


